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I.  Introduction 


The  greatest  advantage  of  a  phased-array  antenna  is  its  capability  to 
scan  a  single  beam  or  multiple  beams  over  a  broad  zingular  region.  The  main 
disadvantage  is  that  to  produce  a  small  beamwidth,  a  large  aperture  and 
thus  a  large  number  of  array  elements  are  required.  The  number  of  elements 
is  determined  by  the  spacing  required  to  avoid  grating  lobes  in  the  field  of 
view.  A  method  used  to  reduce  the  number  of  elements  required  is  to  magnify 
the  size  of  the  array  optically  until  the  beamwidth  requirement  is  satisfied. 
This  is  done  by  using  the  array  as  the  source  feed  to  a  reflector  system  [1-3]. 
Unfortunately,  the  maximum  scan  angle  and  the  grating  lobe  separation  of 
the  combined  system  are  both  inversely  proportional  to  the  magnification. 
Therefore,  there  are  practical  limits  to  the  amount  of  magnification  allowed. 
As  the  array  is  scanned  off  the  boresight  direction,  optical  abberations  occur 
that  distort  the  magnified  image  of  the  array.  The  optical  path  length  and 
the  actual  angle  at  which  radiation  from  the  image  points  leaves  the  system 
are  different  for  each  point.  This  leads  to  phase  error,  which  reduces  the 
peak  gain  of  the  antenna. 

The  scanning  capability  of  an  offset-feed  near-field  Gregorifui  reflector 
is  limited  primarily  by  pha.se  errors  across  the  image  aperture.  A  linear 
phase  across  the  array  plane  produces  a  nonlinear  phase  across  the  image 
plane,  which  increases  with  the  arraj'  scan  angle.  Since  the  reflector  is  fed 
by  a  phased  array,  it  is  possible  to  compensate  for  optical  abberations  by 
adjusting  the  relative  phase  of  the  signal  supplied  to  each  array  element 
|4-7]. 

The  purpose  of  this  study  is  to  demonstrate  a  method  of  correcting  the 
scanning-phase  error  of  a  near-field  offset-feed  Gregorian  reflector  system. 
This  method  views  the  reflector  as  a  simple  magnification  system  that  images 
each  array  element  from  a  point  on  the  array  plane  to  a  i>oint  on  the  leflectoi 
aperture  plane.  Thus,  the  reflector  aperture  plaiie  is  the  same  as  the  image- 
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array  plane.  Requiring  a  linear  phase  across  the  image  array  implies  that 
the  angle  of  each  ray  leaving  the  image  array  aperture  is  the  same.  Thus, 
a  ray  leaving  each  antenna  element  in  the  array  plane  is  traced  through 
the  reflector  to  the  image  aperture-plane  to  find  the  final  exit  angle.  In  a 
perfect  optical  system,  that  angle  should  be  inversely  proportional  to  the 
magnification  of  the  system.  If  it  is  not,  then  the  array-plane  scan  angle  is 
adjusted  imtil  the  exit  angle  matches  the  theoretical  value.  This  produces  a 
scan  phase  for  each  array  element,  together  with  a  phase  that  matches  each 
optical  path  length  to  correct  the  phase  distortion  in  the  image  aperture- 
plane.  Higher-order  effects  such  as  diffraction  and  mutual  coupling  are  not 
considered  for  this  analysis. 
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II.  Principle  of  Reflection 


The  principle  of  reflection  is  somewhat  misleading  because  it  is  based 
on  human  perception  rather  than  on  an  understanding  of  the  physics.  It 
is  “esthetic”  to  view  a  light  ray  as  something  that  bou.ices  from  surface 
to  surface.  However,  since  light  is  part  of  the  electromagnetic  spectrum,  it 
is  well  known  that  this  does  not  adequately  explain  the  phenomenon.  An 
incident  electric  field  (light  ray)  impinging  on  a  surface  induces  volume  and 
surface  currents  to  flow  in  the  object,  which  in  turn  produces  a  secondarj' 
source  of  radiation.  This  process  takes  place  so  fast  that  to  the  human  eye, 
it  appecirs  that  the  light  simply  bounces  from  the  surface. 

Designs  for  reflecting  antennas  aie  derived  from  the  optical  projrerties  of 
parabolic,  hyperbolic,  and  elliptical  reflecting  surfaces.  Using  optical  design 
methods  is  valid  when  the  diameter  of  the  reflector  is  many  times  that  of  the 
wavelength  of  the  radiation.  The  total  electric  field  vector  representing 
the  interaction  of  a  light  ray  with  a  conducting  surface,  can  be  thought  of 
as  the  sum  of  an  incident  vector  field  E'  and  a  reflected  vector  field  such 
that 


E^  =  E•  +  £^  (1) 

The  incident  and  reflected  electric  fields  can  be  decomposed  into  a  tangential 
and  normal  component  with  respect  to  the  reflecting  surface  as 

E'  =  E;  +  E|  (2) 

E^  =  E;  +  A^  (3) 

The  boundary  conditions  on  the  total  electric  field  at  the  surface  of  a  perfect 
conductor  are  given  by 


n  X  E^  =  0 
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(4) 

(5) 


where  p,  is  the  surface  ch2Lrge  on  the  conductor  and  t  is  the  permittivity  of 
the  medium.  Applying  eq.  (4)  to  eq.  (1)  yields  the  relationship  between  the 
tangential  incident  and  reflected  fields  as 

^  =  -E\.  (6) 

This  implies  that  the  incident  and  reflected  tangential  fields  have  equal  mag¬ 
nitudes  but  opposite  signs  at  the  reflector  surface.  The  principle  of  conserva¬ 
tion  of  energy  maj'  be  applied  to  reason  that  the  magnitude  of  the  incident 
field  and  the  magnitude  of  the  reflected  field  must  be  equal  at  the  conduc¬ 
tor  surface.  This  then  implies  that  the  normal  components  must  also  have 
equal  magnitudes  at  the  conductor  surface.  It  is  well-known  that  the  surface 
charge  on  a  perfectly  conducting  surface  illuminated  by  an  incident  wave  is 
nonzero,  and  therefore  the  normal  components  of  the  incident  and  reflected 
fields  must  have  the  same  magnitude  and  direction  at  the  conductor  surface. 
The  relationship  between  the  incident  eind  reflected  electric  fields  may  be 
expressed  in  vector  form  as 


«  X  £•’■  = -(ri  X  F)  (7) 

n  F^ri  F.  (8) 

Taking  the  cross  product  of  eq.  (7)  with  h  and  using  the  vector  identity 

Ax{B  xC)=  B(A  C)-C(A  B)  (9) 


yields  the  expression 

F  =  2(nF)n-F  (10) 

commonly  known  as  Snell’s  law  of  reflection. 

Let  the  propagation  direction  of  the  incident  and  reflected  fields  be  de¬ 
noted  by  the  unit  vectors  s’  and  s’’,  respectively.  Snell’s  law  iinjilies  that  the 
angle  between  the  incident  wave  and  the  surface  normal  is  the  same  as  the 
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angle  between  the  reflected  wave  and  the  normal.  This  implies  that 


(11) 

(12) 


n  X  s’'  =  n  X  s* 
h  ■  s’’  =  — (n  •  s’) 

which  then  leads  to  the  relation 

s’ =s’-2(ns‘)H.  (13) 

This  equation  can  be  put  in  the  form 

3 

—  y  Tmn  (14) 

n  =  l 

where  T^n  represents  the  reflection-matrix  coefficient  taken  from  the  matrix 
equation 

s’  1  —  2n^  — 2Hi.»iy  — 2n,ri.  sj 

s’  =  -2rijny  1  -  2»)'  -2riyTn  sj,  (15) 

2nxn..  — 2nyn;  1  —  2nj  s’. 

Thus,  the  reflection  properties  of  a  perfectly  conducting  surface  can  be  com¬ 
puted  once  the  surface  normals  are  found. 
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III.  Offset-Feed  Near-Field  Gregorian  Design 

The  general  form  of  a  classical  Gregorian  optical  reflector  system  is  shown 
in  Fig.  1.  Note  that  a  spherical  wave  emanating  from  focus  1  of  the  ellipse 
reflects  off  the  elliptical  subreflector  and  is  focused  at  focus  2.  The  spherical 
wave  emanating  from  focus  2  then  reflects  off  the  parabolic  surface  of  the 
main  reflector,  producing  a  plane  wave.  A  typical  antenna  system  would  have 
a  horn  either  receiving  or  radiating  at  focus  1.  The  subreflector  is  assumed 
to  be  in  the  far-zone  field  of  the  radiating  source.  Note  that  the  subreflector 
blocks  a  portion  of  the  reflector  aperture;  this  blockage  is  inherent  in  this 
axially  symmetric  design  and  leads  to  a  loss  in  gain  and  an  increase  in  the 
sidelobe  level. 

If  the  spherical-wave  point  somce  is  replaced  by  a  two-dimensional  array 
of  point  source  radiators,  then  even  though  the  subreflector  is  in  the  far 
zone  of  each  radiating  element,  it  is  in  the  near  zone  of  the  array.  The  wave 
emanating  from  the  array  can  be  appioximated  by  a  plane  wave.  Given 
this  assumption,  the  elliptical  subreflector  must  be  replaced  by  a  parabolic 
subreflector  so  that  the  jdane  wave  produced  by  the  array  can  be  converted 
into  a  spherical  wave  by  the  subreflector.  If  the  array  is  laterally  offset, 
then  it  is  advantageous  to  use  an  offset  geometry  such  as  that  shown  in 
Fig.  2,  to  avoid  any  blockage.  This  offset  configuration  uses  two  confocal 
paraboloidal  reflector  sections.  The  array  configuration  chosen  is  shown  in 
Fig.  3.  The  geometry  is  composed  of  a  hexagonal  lattice  truncated  to  form 
an  approximately  circular  aperture. 
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Fig. 

1.  Gregorian  Dual  Reflector 

Antenna 

Assume  the  array  is 

located  in  the  plane  (x.y,;,,). 

The  eciuation  of  the 

subreflector  is  given  by 

‘  +  S'*) 

(16) 

and  the  equation  of  the  main  reflector  is 

(17) 

Let  R  denote  the  position  vector  on  the  surface  of  either  reflector  given  by 

R  =  jrir  +  yy  +  zt. 

(18) 

/ 

/ 

/ 

/ 

i 

/ 


Fig.  2.  Offset  Gregorian  Dual-Reflector  Antenna 
A  unit  vector  that  is  normal  to  the  surface  may  nov  be  defined  by  the  relation 

( dk  ^ 

\di  ^  dy  ) 


ii  zz  ±. 


dR  ^  <w 

dz  dy 


(19) 


(20) 


which  yields  the  subreflector  unit  normal  n,,  given  by 

-  ~  y.v  - 

+  y=  AF- 

and  the  main  reflector  unit  normal  nm,  given  by 

^  _  -XT  -yy  +  2F„: 

where  the  sign  is  chosen  such  that  the  angle  between  the  normal  and  the 
reflected  ray  is  acute. 


(21) 
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Fig.  3.  Hexagonal  Array  with  C13  Elements 

Assume  that  the  array  is  scanned  only  in  the  t  -  z  plane  such  that  a  scan 
angle  q  can  be  defined  by  the  relation 


taii(o)  =  — — —  (22) 

■Jl  —  ‘0 

where  the  angle  is  positive  when  mecisured  counterclockwise  from  the  positive 
z  axis  in  the  x  -  z  plane.  A  vector  from  the  array  plane  to  the  subreflector 
surface  can  be  defined  as 


^01  =  (xi  -  Xo)r  +  (xi  -  zo)i 

Substituting  eqs.  (IC)  and  (22)  into  eq.  (23)  yields 

1  ' 

Rui  =  z,  -  Za  -  ~{x\ y\)  (lan(u)x  +  ;) 


(23) 


(24) 
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Equating  the  x  components  of  eq.  (24)  yields  the  quadratic  equation  in  xj 


Aix\  +  Bixi  +  Cl  =  0  (25) 

where 

Ai  =  tan(a)  (26) 

Bi  =  4F,  (27) 

Cl  =  -4xoF,  +  {yl  -  4(2,  -  2<.)F,)  tan(a).  (28) 


This  equation  will  have  only  one  feasible  solution  for  xj.  The  coordinates 
of  fho  ray  intercept  on  the  subreflector  aie  then  given  by 


(  X(|  a  =  0 

(29) 

yi  =  yo 

(30) 

=-»-^(a'i  +  y]) 

(31) 

Ei'ch  coordinate  (xi,j/i,2i)  can  then  be  checked  to  make  sure  it  intercepts  the 
physical  subreflector. 

A  general  ray  between  the  subreflector  and  the  main  reflector  can  be 
defined  as 


^12  =  (j‘2  -  a-ilj"  +  (y2  -  yi)y  +  (-2  -  -i)^  (33) 

where  (x2,y2.^2)  is  the  luiknown  point  on  the  main  reflector  surface.  An 
application  of  Snell’s  law  on  the  surface  of  the  subreflector  shows  that  this 
ray  must  be  parallel  to  the  ray  S  defined  by 


5  =  .^1  -  2(n,  •  Roi)n,. 


(34) 


In  order  for  eqs.  (33)  and  (34)  to  be  parallel,  the  slopes  measured  in  two 
orthogonal  planes  must  be  equal.  This  leads  to  the  two  equations 


y2  -  yi  _  ^ 
a’  2  -  ri  Sr 
>2  -  Zi  _  ^ 
-  r,  Sr 


(35) 

(36) 
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where 


Sz  =  T,\\Rxi  +  T,x2Ryi  ATtxzRii 

(37) 

Sy  —  T,2lRzl  +  T,22Ryl  ■AT,23R21 

(38) 

S2  =  T,3i-Ri1  +T,32Ryi  4-7',337?j1. 

(39) 

Combining  eqs.  (35)  and  (36)  yields  a  quadratic  equation  in  X2  given  by 


A2x\  +  B2X2  +  C2  —  0 


(40) 


where 


A2  ~  1 

(41) 

8-2  =  2ui(yj  -  uixi)  -  4u2^m 

(42) 

C2  =  (!/i  -  tiixi)^  +  4F„,(;m  -xi  +U2X1) 

(43) 

s. 

(44) 

5.. 

(45) 

The  coordinates 
given  by 


(i'2.!/j,-:2)  of  the  ray  intercept  on  the  main  reflector  are  then 


^•2 


J-i 


-Bi±y/B^-4A3Cj 

2-42 


Q  =  0,  X()  =  y„  =  0 
a  ^  0 


(46) 


and 


y2  =  yi  +  (a-2  -  a'Oui  (47) 

i2  =  +  yl)  (48) 

The  sign  ambiguity  in  eq.  (46)  is  due  to  the  ambiguity  in  the  ratio  of  the 
slopes.  The  correct  sign  is  chosen  by  making  sure  that 


sign[x2  -  X,]  =  sign[5,] 
sign[y2  -  Vi]  =  sign[5j,]. 
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(49) 

(50) 


In  a  manner  similar  to  that  for  the  subreflector,  each  coordinate 

can  be  checked  to  make  sure  it  intercepts  the  physical  surface  of  the  main- 

reflector. 

The  general  ray  between  the  main  reflector  and  the  aperture  plane  is 
given  by 

R-23  =  {X3  -  X2)x  +  (ys  -  !/2)y  («3  “  ^2)^  (51) 


Applying  Snell’s  law  on  the  surface  of  the  main  reflector  yields 


ys  -  y-2 

—  t’l 

(52) 

23  —  Z2 

X3  -  2-2 

—  f2 

(53) 

ya  -  y-i 

where 

^  Tm3\Rzl  +  Tm^^Ryl  +  TmisRil 

TnillRzi  +  TmliRyi  +  Tmisfizl 

^04/ 

^  7m3ll?rl  4-7m32-Ryl  -l-Tmaa/?.) 

+7'mJ2l?!/l  -1- 7ml3^.-l 

Solving  eqs.  (52)  and  (53)  yields  the  coordinate  point  (x3,y3,23)  of  the  inter- 

ception  of  the  ray  and  the  aperture  plane  of  the  main 

given  by 

reflector;  this  point  is 

^3  =  Xp 

(56) 

ys  =  y2  +  (‘3  “•  •zy)*  1 

(57) 

X3  =  X2  (ys  -  y2)'’2 

(58) 

where  Zp  is  a  specified  aperture  coordinate.  The  total  path  length  of  a  ray  P, 
is  the  total  length  from  the  array  plane  to  the  aperture  plane  and  is  computed 

from  the  formula 

Pi  =  (^.|-H^i2!-t-|^23|. 

(59) 
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IV.  Small  Scan-Angle  Transfer  Function 


The  reflector  scan  angle  0  has  a  complicated  nonlinear  dependence  on 
the  array  scan  angle  a.  However,  a  simple  linear  transfer  function  can  be 
derived  for  small  a  angles.  This  simple  relationship  will  be  used  in  the 
correction  process  to  estimate  the  location  of  the  main  beam.  The  derivation 
is  accomplished  by  taking  a  single  axial  ray,  tracing  it  through  to  the  aperture 
plane,  and  then  making  small  angle  approximations. 

Reflection  from  the  subreflector  is  illustrated  in  Fig.  4.  The  normal 
unit  vector  can  be  found  by  first  noting  that  the  family  of  curves  of  the 
subreflector  can  be  written  as 

J.2 

f(x,  z)  =  z  -  =  positive  constant .  (60) 

4c, 

The  inward  normal  vector  is  defined  by 


n. 


IL 

l^^/l 


which  yields 


n,  = 


(61) 


(C2) 


It  is  useful  to  define  an  angle  0  between  the  unit  vector  h,  and  the  horizontal 
vector  -2  at  the  point  of  reflection  (xi.zj),  which  may  be  represented  as 


=  (63) 

Now  define  (,  as  the  angle  between  the  reflected  ray  and  the  unit  vector  -i 
such  that  ^  =  2/?  +  Q. 
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surface  normal 


Fig.  4.  Reflection  from  Subreflector 

Figure  5  illustrates  reflection  from  the  main  reflector  surface.  The  coor¬ 
dinates  of  a  ray  interception  on  the  main  reflector  surface  are  given  by 


tan^ 


T2  —  Ti 
Zl  -  22 


which  leads  to  the  quadratic  equation 


(64) 


4”  ^'■1^2  4"  6^2  —  0 


(65) 


where 


A2  = 


4Fm 


(66) 


B2  =  1  (67) 

C2  =  (F.  tan(0  -  *i  (68) 


The  coordinates  of  the  ray  intercept  on  the  main  reflector  are  given  by 


T2 


Xi 

—  fij—  y/  — 4^  jCj 

2/t, 
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f  =  0 


(69) 


A  unit  vector  normal  to  the  main  reflector  pointing  in  the  direction  of 
the  subreflector  is  given  by 


n  - 

»‘m  —  « 


Then  defining  an  angle  y  between  the  imit  vector  i  and  at  the  reflection 
point  (x2,22)  yields  the  relation 


7  =  tan 


Defining  an  angle  0  as  the  angle  that  the  ray  reflected  off  the  main  reflector 
makes  with  the  unit  vector  2,  such  that  5  =  4  -  27,  it  follows  that 


tan(d) = 


X3  - 
23  -  22 


which  can  be  rewritten  in  terms  of  the  aperture-plane  t  coordinate  as 


23  =  2'2  +  (-3  -  22)tan(/?) 


The  reflector  scan  angle  6  and  the  array  scan  angle  o  are  related  as 


6  =  a  +  2{l3-y). 


For  small  scan  angles  this  relationship  may  be  approximated  as 


0  «  /(O)  -I-  /'(0)a 


where 


/(a)  =  Q  +  2[d(<»)  -  7(0)] 


Fig.  5.  Reflection  fiom  Main  Reflector 
By  observation  it  can  be  seen  that  /(O)  =  0  such  that 

6  «  /'(0)q 

where 


da 


Using  eq.  (71),  it  can  be  shown  that 


df 


2F„ 


da  4F^  + 

Using  implicit  differentiation  on  eq.  (65)  yields 


j  Ml  .,2  _i_  ic, 

^  _  ~dS^^2~*'  ~d^ 

da  2A2T2  "f  ^2 


where 


dAo 


l  +  2f. 


da  cos2(o  +  20) 
da  V'  4fJ 

Evaluated  at  a  =  0  yields 


x  +  24^ 


da 


cos^(2;3 


_  1  ^ 

2  <fo 


iL_^_fl  +  ^tan(2/?  +  a)')$l 
+  a)  V  2F.  '  7  do 


171 

2F„ 


(77) 

(78) 

(79) 

(80) 

(81) 

(82) 

(83) 
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Substituting  this  into  eq.  (78)  yields 

/'(0)  =  -p-  (84) 

Therefore,  for  small  array  scan-angles,  the  array  scan  angle  a  and  the  reflector 
scan  angle  are  related  as 

e  =  -^o.  (85) 

*  m 

This  indicates  that  the  maximum  reflector  scan  angle  is  inversely  propor¬ 
tional  to  the  magnification  of  the  system.  This  implies  that  the  scanning 
capability  of  a  reflector  system  is  directly  related  to  the  magnification. 
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V.  Small  Scan-Angle  Ray-Bundle  Characteristics 

The  side  view  of  the  near-field  Gregorian  system  with  a  boresight  beam 
is  shown  in  Fig.  6.  The  array  has  a  diameter  of  Da  =  ISAq.  A  hexagonal 
array  of  613  elements  fits  within  the  diameter,  with  a  center-to-center  element 
spacing  of  s  =  .6923Ao.  The  array  is  positioned  in  the  x~y  plane  at  i  =  0  with 
a  horizontal  offset  position  of  Xa  =  -9Ao-  The  subreflector  and  main  reflector 
have  focal  lengths  of  F,  =  24Ao  and  =  96 Ao,  for  a  magnification  oi  M  =  A. 
The  on-axis  offsets  are  given  by  z,  =  42Ao  and  z^  =  -78Au.  The  subreflector 
offset  limits  are  given  by  xt  =  -40Ao  and  x^  =  6.2Ao.  The  main-reflector  offset 
limits  are  given  by  xj  =  15Ao  and  x^  =  llOAo.  The  rays  are  terminated  at  the 
reflector  aperture  plane  located  at  Zp  =  -46  49.  For  this  scan  angle,  all  rays 
pass  through  the  focus  and  form  an  exact  point  source.  The  front  vdew  of  the 
ray  bundle  is  shown  in  Fig.  7.  Each  array  point  of  the  source  is  imiformly 
magnified  and  reflected  across  the  center  of  the  array.  There  is  no  optical 
distortion. 

The  side  view  of  the  uncorrected  ray  bundle  for  a  =  -27®  is  shown  in  Fig. 

8.  The  focus  point  moves  below  the  axis  and  is  blurred.  Each  ray  comes 
out  of  the  main-reflector  aperture  plane  at  a  different  angle  but  near  the 
design  angle  o{  0  =  6.75°.  The  front  view  of  the  ray  bundle  is  shown  in  Fig. 

9.  Note  that  the  equivalent  area  of  the  array  is  compressed  in  the  vertical 
direction.  This  will  increase  the  beamwddtli  slightly  in  the  vertical  plane. 
For  the  applications  of  interest  in  this  study,  the  increased  beamwidth  is  not 
considered  critical. 

The  side  view  of  the  uncorrected  ray  bundle  for  a  =  20°  is  shown  in  Fig. 

10.  The  focus  point  moves  above  the  axis  and  is  blurred.  Each  ray  comes  out 
of  the  main  reflector  aperture  plane  at  a  different  angle  but  near  the  design 
angle  of  0  =  -5°.  The  front  view  of  the  ray  bundle  is  shown  in  Fig.  11.  Note 
that  the  equivalent  area  of  the  array  is  expanded  in  the  vertical  direction. 
This  will  decrease  the  beamwidth  slightly  in  the  vertical  plane. 
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Fig.  10.  Side  View  of  Raj'  Bundle  with  a  =  20° 
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VI.  Scanning-Phase  Correction 


The  two  main  types  of  phase  errors  are  caused  by  unequal  path  lengths 
and  unequal  scan  angles  for  rays  that  originate  in  the  array  plane  and  ter¬ 
minate  in  the  aperture  plane.  Each  of  these  errors  can  be  eliminated  by 
appropriately  adjusting  the  phase  on  each  array  element.  At  each  scan  an¬ 
gle,  a  set  of  correction  phases  can  be  computed  £ind  used  to  cancel  any  phase 
distortion. 

Because  of  the  nonlinear  relationship  between  a  and  a  uniform  scan 
angle  applied  to  the  array  plane  yields  a  nonuniform  scan  angle  across  the 
main-reflector  aperture  plane.  This  gives  a  pointing-beam  error  as  well  as  a 
loss  in  gain.  It  is  well-known  that  a  linear  phase  is  required  to  scan  a  beam 
with  no  loss  in  power.  Therefore,  each  raj-  must  leave  the  main  reflector  at 
the  same  angle  0.  The  value  of  0  is  derived  from  the  small  scan-angle  relation 
given  in  eq.  (85).  Using  this  as  a  target,  the  nth-element  array  scan  angle  q„ 
lies  within  the  range 

^0-Aa<a„<  +  (86) 

r,  F, 

A  simple  bisection  search  algorithm  is  performed  over  the  interval  until  the 
element  scan  angle  yields  the  desired  reflector  scan  angle  0. 

Once  the  scan  angle  correction  has  been  made,  the  path  length  from  the 
array  plane  to  the  aperture  plane  can  be  computed.  The  path  length  for  a 
scan  angle  of  q  =  0°  is  constant  and  given  by 


Fo  —  2(*.  Za)  "t"  (Za 


(87) 


As  the  scan  angle  changes,  the  path  length  for  each  ray  changes.  The  path 
length  error  in  radians  is  given  by 


A,  =  2w 


(Pi  -  Fo) 
^0 


(88) 


This  phase  correction  is  added  to  each  array  element,  so  that  each  ray  appears 
to  have  traveled  the  same  distance. 


33 


VII.  Results 


Since  the  antenna  system  is  designed  so  that  there  is  no  spillover  for 
either  reflector,  the  correction  phase  for  each  array-element  ray  should  be 
constrained  to  allow  no  spillover.  Fortunately,  it  turns  out  that  the  maximum 
correction  angle  for  all  rays  for  negative  a  is  greater  than  -27®.  For  positive  q 
the  maximum  correction  angle  is  several  degrees  greater  than  20®.  However, 
the  rays  that  need  the  most  correction  are  those  with  the  greatest  horizontal 
offset,  so  even  though  the  corrected  angle  is  greater  than  20°,  the  ray  still 
intercepts  the  subreflector  with  no  spillover. 

The  relationship  of  the  design  arraj'  scan  angle  to  the  reflector  scan 
angle  for  the  uncorrected  and  corrected  cases  is  shown  in  Fig.  12.  The 
reflector  scan  angle  is  defined  to  be  the  angular  point  of  the  main-beam 
peak.  The  corrected-case  line  has  a  slope  of  -4,  equal  to  the  negative  of  the 
magnification.  The  nonlinear-  imcorrected  curve  shows  the  increase  in  error 
as  the  array  scan  angle  moves  from  zero.  The  main-beam  peak  power  as  a 
function  of  the  array  scan  angle  is  shown  in  Fig.  13.  Note  that  the  corrected 
curve  indicates  no  loss  in  main-beam  power  over  the  entire  scan  range. 

These  two  figures  indicate  that  the  main-beam  peak  power  is  maximum 
and  constant  when  the  reflector  scan  angle  is  linearly  related  to  the  design 
array  scan  angle.  Full  correction  can  be  acliieved  by  forcing  the  reflector 
scan  angle  to  be  uniform.  However,  it  should  be  noted  that  changes  to 
the  beamwidth  and  grating  lobe  separation  axe  not  restricted.  The  change 
in  beamwidth  and  grating  lobe  separation  could  be  minimized  by  using  a 
shaped  subreflector. 
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VIII.  Conclusion 


Scanning-phase  error  for  an  offset-feed  near-field  Gregorian  reflector  can 
be  corrected  by  adjusting  the  phase  of  each  array  element.  A  nonlinear 
array  scan-angle  distribution  can  be  found  that  yields  a  hnear  reflector  scan 
angle  that  is  the  same  for  each  element  beam.  This  phase  correction  cam  be 
found  by  a  numericail  search  in  the  vicinity  of  the  linear  design  amgle.  For  the 
geometrical  parameters  considered,  the  phase  correction  can  be  accomplished 
without  any  loss  in  main-beam  power  or  scan  range. 
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Appendix:  Computer  Program  PHGREG 


This  appendix  contains  the  Fortran  77  program  PHGREG,  tis  well  as  all 
the  associated  subroutines.  The  program  implements  the  scanning  phase- 
correction  procedure  described  in  this  report.  All  input  and  output  variables 
are  defined  in  the  comment  section  at  the  beginning  of  the  program  and  each 
subroutine. 
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PROGRAM  PHGREG 

*******t********************************************************** 


1 

2  C 

3  C 

4  C 

5  C 

6  C 

7  C 

8  C 

9  C 

10  C 

11  c 

12  C 

13  C 

14  C 

15  C 

16  C 

17  C 

18  C 

19  C 

20  C 

21  C 

22  C 

23  C 

24  C 

25  C 

26  C 

27  C 

28  C 

29  C 

30  C 

31  C 

32  C 

33  C 

34  C 

35  C 

36  C 

37  C 

38  C 

39  C 

40  C 

41  C 

42  C 

43  C 

44  C 

45  C 

46  C 

47  C 

48 

49 


THIS  PROGRAM  SYNTHESIZES  THE  SCANNING  PHASE  OF  AN  OFFSET-FEED  * 
NEAR-FIELO  GREGORIAN  REFLECTOR.  * 

**«******4*********««****4*4«*********4*********4«**«****4******* 

TIMOTHY  J.  PETERS  LAST  UPDATED  ♦ 

THE  AEROSPACE  CORPORATION  1/15/92  * 

2350  EAST  EL  SEGUNDO  BOULEVARD.  * 

EL  SEGUNDO,  CA  90245  * 

**m*t***t********************«*********************************** 

INPUTS:  • 

* 

MN  -  MAXIMUM  ARRAY  DIMENSION.  * 

DA  -  ARRAY  DIAMETER  IN  WAVELENGTHS.  * 

ZA  -  ALONG  AXIS  ARRAY  LOCATION  IN  WAVELENGTHS.  * 

XA  -  LATERAL  DISPLACEMENT  OF  ARRAY  IN  WAVELENGTHS.  * 

ZS  -  ALONG  AXIS  SUBREFLECTOR  LOCATION  IN  WAVELENGTHS.  ♦ 

XB  -  MINIMUM  LATERAL  OFFSET  OF  SUBREFLECTOR  GEOMETRY  IN  • 

WAVELENGTHS .  * 

XC  -  MAXIMUM  LATERAL  OFFSET  OF  SUBREFLECTOR  GEOMETRY  IN  * 

WAVELENGTHS .  * 

FS  -  FOCAL  LENGTH  OF  SUBREFLECTOR  IN  WAVELENGTHS.  * 

FM  -  FOCAL  LENGTH  OF  MAIN  REFLECTOR  IN  WAVELENGTHS.  * 

XD  -  MINIMUM  LATERAL  OFFSET  OF  MAIN  REFLECTOR  GEOMETRY  IN  • 

WAVELENGTHS .  * 

XE  -  MAXIMUM  LATERAL  OFFSET  OF  MAIN  REFLECTOR  GEOMETRY  IN  • 
WAVELENGTHS .  • 

ZM  -  ALONG  AXIS  MAIN  REFLECTOR  LOCATION  IN  WAVELENGTHS.  ♦ 

NL  -  NUMBER  OF  COMPLETE  HEXAGONAL  FOR  ARRAY  GEOMETRY  ♦ 

GENERATION .  * 

ISIN  -  ISIN=0  MEANS  NO  CORRECTION  IS  PERFORMED  AND  ISIN=1  ♦ 

MEANS  FULL  CORRECTION  OF  PHASE  ERROR.  * 

SMIN  -  MINIMUM  AND  MAXIMU  ARRAY  SCAN  ANGLES  IN  DEGREES.  ♦ 

SMAX  * 

NSA  -  NUMBER  OF  SCAN  ANGLES.  ♦ 

TOL  -  ERROR  TOLERANCE  FOR  BISECTION  ALGORITHM.  * 

MI  -  MAXIMUM  NUMBER  OF  ITERATIONS  FOR  BISECTION  ALGORITHM.  * 

* 

OUTPUTS:  * 

4 

SCMIN  -  MINIMIM  AND  MAXIMUM  CORRECTION  ANGLES  IN  DEGREESFOR  ♦ 
SCMAX  FOR  A  GIVEN  ARRAY  SCAN  ANCLE  WRITTEN  TO  FILE  'RANGE’.* 
RPEAK  -  MAIN  BEAM  PEAK  ANGLE  IN  DEGREES  FOR  A  GIVEN  ARRAY  * 
SCAN  ANGLE  WRITTEN  TO  FILE  'ANGLE'.  * 

PDB  -  MAIN  BEAM  PEAK  POWER  IN  DB  FOR  A  GIVEN  ARRAY  SCAN  * 
ANGLE  WRITTEN  TO  FILE  'GAIN'.  * 


**********************«******************************4******444*** 

PARAMETER  (MN=1000) 

REAL*4  X(MN) .Y(MN) ,XT(MN) ,YT(MN) .ASA(MN) . V(MN) 
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so  IHTEGER*4  IRC(MN) 

51  PI=3.141S93 

52  RAD= . 17453293E-01 

53  OPEH(UMIT= 10, F1LE=’ RANGE *) 

54  OPEN (U«IT=11.FILE=* ANGLE*) 

55  0PEN(UNIT=12,FILE=’GAIN’) 

56  C  «4««*44*******44*«««4*4«««*««««***»**»**44*444****4*«4*44**4**«**4 

57  C  ♦  REFLECTOR  GEOMETRY  INPUTS.  * 

58  C  *****«*4*******«*«**«***«**«*******«****»4«************4>*»4*«***** 

59  DA=18.0 

60  ZA=0.0 

61  XA=-9.0 

62  ZS=42.0 

63  XB=-40.0 

64  XC=6.2 

65  FS=24.0 

66  FM=96.0 

67  XD=1S.0 

68  XE=110.0 

69  ZM=-FM+2S-FS 

70  NL=13 

71  ISIN=1 

72  SMIH=-27.0 

73  SMAX=20.0 

74  NSA=95 

75  T0L=O. 000001 

76  HI=80 

77  C  «**8**4>44**4*****4'*4'*****4************4*4i**4*****44**44***4******* 

78  C  *  GENERATE  THE  ARRAY  ELEMENT  GEOMETRY.  ♦ 

79  C  ****♦**♦♦♦♦♦**♦** ♦♦♦♦♦♦*4************4*************»4*****»*****»4 

80  RA=DA/2.0 

81  S=RA/NL 

82  XAC=XA-DA/2.0 

83  YAC=0.0 

84  CALL  GEOMET(MN,NL.S.XAC,YAC.X,Y,NE) 

85  C  ♦♦♦♦♦**♦**♦*•♦♦♦♦♦♦♦♦♦*♦♦**♦****♦•••**♦♦♦♦♦♦*♦♦♦****»♦♦♦**»****♦♦♦ 

86  C  *  COMPUTE  THE  PEAK  POWER  POINT  AT  BORESIGHT.  • 

87  C  4 4 4444* 4444**4*4 ****4*«4***4******** ******* 4 *444 *44 4* 4* ******** 

88  DO  1  I-1,NE 

89  ASA(I)=0.0 

90  V(I)=0.0 

91  1  CONTINUE 

92  CALL  PEKPOW(MN,NE,X.Y.XT,YT.ASA.V.ZA,FS.ZS,FM,ZM.XB,XC.XD.XE 

93  ft  ,IRC,ISIN.RPEAK,PMAXG) 

94  C  44****44444*444**4*44*44***4*44*4*****444444*4*4**44***4444**4**4* 

95  C  ♦  STEP  THROUGH  THE  SCAN  ANGLES.  * 

96  C  44444*44*444**4*44*44444*4*4*4*444*****44444*44**444*44*44*4*44444 

97  DSA=(SMAX-SMIN)/(NSA-1) 

98  DO  2  K=0,NSA-1 
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99 

100  C 

101  C 

102  C 

103 

104 

105 

106  C 

107  C 

106  C 

109 

110 

111  C 

112  C 

113  C 

114 

115 

116 
117  C 
116  C 

119  C 

120  C 

121 
122 

123 

124 

125 

126 
127 
126 

129 

130 

131 

132 

133 

134 

135 

136 

137 
136 

139 

140  4 

141  99 

142  C 

143  C 

144  C 

145 

146 

147  3 


SCAK=SMIH+K*DSA 

66*444*6***** *6 *4**«**«**»*«****6*******46**6«*6**4>«44*****4**** 

*  GERERATE  THE  SCANNIKG  PHASE  APPLIED  TO  EACH  ARRAY  ELEMENT.  * 

**************************************************************** 

TARGET=-(FS/FM)*SCAN 
DO  3  1=1, RE 

IF  (ISIN  .EQ.  0)  THEN 

************************************************************ 

*  RO  CORRECilOR.  • 

************************************************************ 

ASA(I)=SCAN 

ELSE 

************************************************************ 

*  WITH  CORRECTION.  * 

************************************************************ 

XO=X(I) 

YO=Y(I) 

ZO=ZA 

************************************************************ 

*  USE  A  BISECTION  SEARCH  ALGORITHM  WITH  +-10  DEGREES  FROM  * 

*  THE  TARGET  SCAN  ANGLE.  • 

************************************************************ 

AMIN=SCAN-10.0 
AMAX=SCAN+10.0 
DO  4  J=1,MI 

A=AMIN+(AMAX-AMIN)/2.0 

CALL  RAYATS(XO,YO,ZO,A,FS,ZS,FM.ZH.XE,X1,Y1,Z1 
t  ,X2,Y2,Z2.X3.Y3,Z3) 

T=X3-X2 
B=Z3-Z2 

THETA=180*ATAN(T/B)/PI 

ERR=THETA-TARGET 

IF  ((ABS(ERR)  .LT.  TOL)  .OR.  (J  .EQ.  MI))  THEN 
GO  TO  99 
ELSE 

IF  (ERR  .GT.  0.0)  THEN 
ANIN=A 
ELSE 
AMAX=A 
END  IF 
END  IF 
CONTINUE 
CONTINUE 

************************************************************ 

*  STORE  THE  COMPUTED  SCAN  ANGLE.  * 

************************************************************ 

ASA(I)=A 
END  IF 
CONTINUE 
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#**4***4**4*4***«**44*««**********4444***4*44***«44************* 

*  FIKD  THE  NIMIMUM  AND  MAXIMUM  ARRAY  SCAN  ANGLES  REQUIRED  TO  * 

*  CORRECT  THE  PHASE  ERROR.  * 

^**t**t********************************************************* 

DO  S  1=1, IE 

IF  (I  .EQ.  1)  THEN 
SCMIN=ASA(I) 

SCMAX=ASA(I} 

ELSE  IF  (ASA(I)  .LT.  SCMIN)  THEN 
SCMIN=ASA(I) 

ELSE  IF  (ASA(I)  .GT.  SCMAX)  THEN 
SCMAX=ASA(I) 

ELSE 
END  IF 
CONTINUE 

***4******44**»44******4***4*******4>4******************4*4**44*4 

*  FIND  THE  PEAK  POWER  LOCATION  AND  VALUE.  * 

«4***«********44*****«***4*«*************»*4***44***4»**»444444* 

CALL  PEKPOW(MN,NE.r,Y,XT,YT,ASA.V,ZA,FS.ZS,FM.ZM.XB,XC,XD,XE 
A  .IRC.ISIN.RPEAK.PMAX) 

**4**«***4**«****4*******4***44***44**4****4*****4****444**4***4 

*  WRITE  OUT  THE  MIN  AND  MAX  SCAN  ANGLES  REQUIRED  FOR  EACH  * 

*  DESIGN  SCAN  ANGLE.  * 

******** **********44«***«»** ******4***** ****»»•*» **•**«******«** 

WRITEClO,*)  K. SCAN, SCMIN, SCMAX 

**************************************************************** 

*  WRITE  OUT  THE  ANGLE  OF  THE  MAIN  BEAM  PEAK  AT  EACH  SCAN  ANGLE.* 

**************************************************************** 

RPEAK=RPEAK/RAD 
WRITECll,*)  SCAN,RPEAK 

**************************************************************** 

*  WRITE  OUT  THE  MAIM  BEAM  POWER  IN  DB  AT  EACH  SCAN  ANGLE.  * 

**************************************************************** 
RAT=PMAX/PMAXG 

PDB=10.0*AL0G10(RAT) 

WRITE(12,*)  SCAN.PDB 
CONTINUE 
CLOSE(IO) 

CLOSE(ll) 

CL0SE(12) 

END 


1  SUBROUTINE  GEOMET(MN.NL,S,XC.yc.X,Y,NE) 

2  C  ***************************************************************** 

3  C  *  THIS  SUBROUTINE  GENERATES  A  CIRCULAR  ARRAY  OF  ELEMENT 

4  C  *  LOCATIONS  BY  TRUNCATING  AN  HEXAGONAL  ARRAY. 

5  C  ***************************************************************** 

6  C  *  INPUTS: 

7  C  * 

8  C  *  MN  -  MAXIMUM  ARRAY  DIMENSION. 
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9  C 
10  C 
n  C 

12  C 

13  C 

14  C 

15  C 

16  C 

17  C 

18  C 

19 
2(1 
21 
22 

23  C 

24  C 

25  C 

26 

27 

28 

29  C 

30  C 

31  C 

32 

33 

34 

35 

36  C 

37  C 

38  C 

39 

40  C 

41  C 

42  C 

43 

44 

45 

46 

47 

48 

49  C 

50  C 

51  C 

52 

53 

54 

55 

56 

57 


♦  NL  -  HUMBER  OF  COMPLETE  UYERS  OF  ELEMENTS.  ♦ 

*  S  -  CENTER  TO  CENTER  ELEMENT  SPACING.  * 

•  XC.YC  -  CENTER  OF  ARRAY.  * 

*  * 

♦  OUTPUTS :  ♦ 

*  * 

♦  X(MN)  -  COORDINATES  OF  THE  CENTER  OF  EACH  ELEMENT.  * 

♦  Y(MN)  ♦ 

♦  NE  -  TOTAL  NUMBER  OF  ELEMENTS.  ♦ 


***«*«*««********«*********«**************48******** 8* ************ 

REAL*4  X(MN).Y(MN) 

REAL*8  RA,W.DP.PHI.CP.SP.T,B.SB.CB.ANGB.CX,SX.XX.YY 

RAD=. 17453293E-01 

PI=.3141593E+01 

****************************************************************** 

♦  SET  THE  CENTER  ELEMENT.  * 

****************************************************************** 

X(l)=XC+0.0 

Y(l)=YC+0.0 

K=1 

****************************************************************** 

♦  LOOP  THROUGH  EACH  LAYER.  * 

****************************************************************** 

RA=NL*S 
DO  1  I=l.NL+5 
W=I*S 

DP=PI/(3*I) 

**************************************************************** 
♦  LOOP  THROUGH  EACH  OF  THE  6  QUADRANTS.  * 

**************************************************************** 

DO  2  L=0,5 

************************************************************** 

•  LOOP  THROUGH  EACH  POINT.  * 

************************************************************** 

PHI=L*RAD*60.0 

CP=DCOS(PHI) 

SP=DSIN(PHI) 

DO  3  J=0,I-1 
T=J*S 

R=SqRT(W*W+T*T-W*T) 

************************************************************ 

*  COMPUTE  ANGLES  FOR  A  PURE  HEXAGONAL  ARRAY.  * 

************************************************************ 

SB= (DSQRT(3 . ODO) /2 . 0)*T/R 
CB=DSQRT(1.0-SB*SB) 

ANGB=DASIN((DSQRT(3.0DO)/2.0)*T/R) 

CX=DCOS(ANGB) 

SX=DSIN{ANGB) 

XX=R*CB 
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56 

59  C 

60  C 

61  C 

62 

63 

64 

65 

66 
67 

66  3 

69  2 

70  1 

71 

72 

73 


YY=R*SB 

«6********6*«**«*****»*******************4*44*464**6***«**6* 

♦  ROTATE  THE  POIMT.  ♦ 

*****««66*****************«****«***«4>****««**************6** 

IF  (R  -LE.  RA)  THEM 
K=K+1 

X(K)=XC+CP*XX-SP*YY 
Y(K)=YC+SP*XX+CP*YY 
ELSE 
EMD  IF 
CONTIMUE 
CONTIMUE 
COHTIHUE 
ME=K 
RETURN 
END 


1  SUBROUTINE  PATHLN(X0.Y0.Z0,X1.Y1.Z1,X2,Y2,Z2.X3,Y3,23,PL) 

2  C  **6«***«*44*********4****6*****************»********»»*****«i»«**^* 

3  C  ♦  THIS  SUBROUTINE  COMPUTES  THE  PATH  LENGTH  OF  A  RAY.  ♦ 

4  C  «*•***«•** •*******«4**********«**«*******»»»»**«*«****»*********»^ 

5  C  *  INPUTS:  * 

6  C  *  • 

7  C  *  (XO.YO.ZO)  -  COORDINATES  OF  POINT  ON  ARRAY.  * 

6  C  *  (Xl.Yl.Zl)  -  COORDINATES  OF  POINT  ON  SUBREFLECTOR.  • 

9  C  *  (X2.Y2.Z2)  -  COORDINATES  OF  POINT  ON  MAIN  REFLECTOR.  * 

10  C  *  (X3,Y3.Z3)  -  COORDINATES  OF  POINT  IN  APERTURE  PLANE.  * 

11  C  *  * 

12  C  *  OUTPUT:  ♦ 

13  C  *  * 

14  C  *  PL  -  PATH  LENGTH.  * 

15  C  ♦♦*♦♦♦♦♦*♦♦*♦♦**•♦**♦»**♦♦*♦*»♦*♦•**♦♦**»*♦**»♦*♦♦♦♦*»♦♦♦♦♦♦**♦♦♦* 

16  DXX=X1-X0 

17  DYY=Y1-Y0 

16  DZZ=Z1-Z0 

19  PLl=SqRT(DXX*DXX+DYY*DYY+DZZ*DZZ) 

20  DXX=X2-X1 

21  DYY=Y2-Y1 

22  DZZ=Z2-Z1 

23  PL2=SqRT(DXX»DXX+DYY*DYY+DZZ*DZZ) 

24  DXX=X3-X2 

25  DYY=Y3-Y2 

26  DZZ=Z3-Z2 

27  PL3=SqRT{DXX*DXX+DYY*DYY+DZZ*DZZ) 

26  PL=PL1+PL2+PL3 

29  RETURN 

30  END 


1  SUBROUTINE  RAYATS(X0,Y0.Z0,AG,FS,ZS.FM,ZM,XE,X1 .Yl.Zl 

2  k  ,X2,Y2,Z2,X3,Y3.Z3) 
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3  C 

4  C 

5  C 

6  C 

7  C 

8  C 

9  C 

10  C 

11  C 

12  C 

13  C 
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15  C 

16  C 

17  C 
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47 

48  C 
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50  C 
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«**«4************************4***«****4*4*4****4****4************* 

♦  THIS  SUBROUTIltE  COMPUTES  THE  INTERSECTION  OF  A  RAY  ORIGINATING  ♦ 


*  ON  THE  ARRAY  VITH  THE  SUBREFLECTOR,  MAIN  REFLECTOR  AND  THE  * 

*  MAIN  REFLECTOR  APERTURE  PLANE.  * 

*****************««*«**«***«************«*****«*«**444**4****«***« 

*  INPUTS ;  ♦ 

«  * 

*  XO.YO.ZO  -  COORDINATES  OF  POINT  ON  ARRAY.  ♦ 

*  AG  -  SCAN  ANGLE  IN  DEGREES.  * 

*  FS  -  FOCAL  LENGTH  OF  SUBREFLECTOR.  ♦ 

*  ZS  -  ALONG  AXIS  POSITION  OF  THE  SUBREFLECTOR.  * 

*  FM  -  FOCAL  LENGTH  OF  MAIN  REFLECTOR.  * 

*  ZM  -  ALONG  AXIS  POSITION  OF  THE  MAIN  REFLECTOR.  * 

*  XE  -  MAXIMUM  LATERAL  EXTENT  OF  MAIN  REFLECTOR.  * 

*  » 

♦  OUTPUTS :  ♦ 

*  * 

♦  Xl.Yl.Zl  -  COORDINATES  OF  POINT  ON  SUBREFLECTOR.  ♦ 

♦  X2.Y2,Z2  -  COORDINATES  OF  POINT  ON  MAIN  REFLECTOR.  ♦ 

♦  X3.Y3,Z3  -  COORDINATES  OF  POINT  IN  APERTURE  PLANE  OF  MAIN  * 

♦  REFLECTOR .  * 

«******4**4**«**«*****«**4****************4***********4******«**** 

RAD=.17463293E-01 

•«««»*«**«**4«********4***«***4**********4*»4»**»*»**»**»*****»»** 

♦  SET  Y1=Y0  SINCE  THE  ARRAY  IS  SCANNED  IN  X-Z  PLANE  ONLY.  * 

**4****4**********************«*************»****»*»*****»4******* 

Y1=Y0 

k***************************************************************** 

♦  COMPUTE  XI.  * 


****4**444*4*444****4****«****4**4*4*************4**********44**** 

ARG=RAD4AG 

TA=TAN(ARG) 

A1=TA 

B1=4.0*FS 

C1=-4.0*X0»FS+(Y1*Y1-4.0*(ZS-Z0)*FS)*TA 
IF  (ABS(Al)  .GT.  0.0)  THEN 
ARG=B1*B1-4.0*A1*C1 
Xl=(-Bl+SqRT(ARG))/(2*Al) 

ELSE 
X1=X0 
END  IF 

*****4**4**4*************4**********4****4444«4*4*********444*4*** 

♦  COMPUTE  Zl.  ♦ 

*4444********4*44****«*44********44**44*********4*****4*4******4*« 

Z1=ZS-(X1*X1+Y1*Y1)/(4*FS) 

*«444*4***4«4«*4***«******4***4*4*44*4*4****4*«***44**4*****4444** 

*  COMPUTE  REFLECTION  MATRIX  FOR  THE  SUBREFLECTOR.  ♦ 

444*** *4 ***4*** ********** *****44**44*4** 4* ************ *4 ********** 

R01X=(Z1-Z0)*TA 


52  R01Y=0.0 

53  R012=Z1-Z0 

54  DEli0H=SQRT(Xl*Xl+Yl*Yl+4.0*FS*FS) 

55  VNX=-X1/DEN0M 

56  VIIY=-Y1/DEH0M 

57  VBZ=-2.0*FS/DE1IQM 

58  TS11=1.0-2.0*VKX*VNX 

59  TS12=-2.0*VHX*VHY 

60  TS13=-2.0*VKX4VKZ 

61  TS21=TS12 

62  TS22=1.0-2.0*VIIY*VHY 

63  TS23=-2.0*VlfY*VlIZ 

64  TS31=TS13 

65  TS32=TS23 

66  TS33=1.0-2.04VNZ*VNZ 

67  C  «4**444«***««***«*******4********«********^********»»**«******««** 

68  C  *  COMPUTE  THE  COOEFICIEHTS  OF  THE  QUADRATIC.  * 

69  C  *44*4*4** *****»44«** ***«***«««**»***#*«*«» *******»4i^»**««« **««***« 

70  XC0MP=TS11*R01X+TS12*R01Y+TS13*R01Z 

71  YCOMP=TS214R01X+TS22*R01Y+TS23*R01Z 

72  ZC0MP=TS31*R01X+TS32*R01Y+TS33*R01Z 

73  C  ****************************************************************** 

74  C  ♦  COMPUTE  THE  POINT  WHERE  THE  RAY  REFLECTED  OFF  THE  SUBREFLECTOR  * 

75  C  *  INTERCEPTS  THE  MAIN  REFLECTOR.  * 

76  C  ****************************************************************** 

77  U1=YC0HP/XC0MP 

78  U2=ZC0MP/XC0MP 

79  A2=1.0+U14U1 

80  B2=2.04Ul*(yi-U:*Xl)-4.0*U2*FM 

81  C2=(Y1-U1*X1)*(Y1-U1*X1)+4.0*FM*(ZM-Z1+U2*X1) 

82  IF  ((ABS(ALPHA) .EQ.  0.0)  . AND . (ABS(XO) . EQ .  0.0) 

83  t  .AND. (ABS(YO)  .EQ.  0.0))  THEN 

84  X2=X1 

85  ELSE 

86  ARG=B2*B2-4*A2*C2 

87  X2A=(-B2-SQRT(ARG))/(2*A2) 

88  X2B=(-B2+SqRT(ARG))/(2*A2) 

89  Y2A=Y1+(X2A-X1)*U1 

90  IF  (Y2A-Y1  .GT.  0.0)  THEN 

91  Pl=1.0 

92  ELSE 

93  Pl=-1.0 

94  END  IF 

95  IF  (YCOMP  .GT.  0.0)  THEN 

96  P2=1.0 

97  ELSE 

98  P2  =  -1.0 

99  END  IF 

100  IF  (X2A-X1  .GT.  0.0)  THEN 
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101 

102 

103 

104 

105 

106 
107 
lOS 

109 

110 
111 
112 

113 

114 

115 

116 
117 

llH  C 

119  C 

120  C 

121 
122 

123 

124 

125 

126 
127 
12» 

129 

130 

131 

132 

133 

134 

135 

136 

137  C 
13S  C 

139  C 

140  C 

141 

142 

143 

144 

145 

146 

147 
146 


P3=1.0 
ELSE 
P3=-1.0 
EHD  IF 

IF  (XCOMP  .GT.  0.0)  THEH 
P4=1.0 
ELSE 
P4=-1.0 
EMD  IF 

IF  ((PI  .EQ.  P2).AND.(P3  .EQ.  P4))  THEN 
X2=X2A 
ELSE 
X2=X2B 
END  IF 
EHD  IF 

Y2=Y1+(X2-X1)*U1 

22=(X24X2+Y2*Y2)/(4«FM)+ZM 

*4*«4**4*4*«*****4********4*******4*4*4**»*****444****«***4**444** 

*  COMPUTE  REFLECTION  MATRIX  FOR  THE  MAIN  REFLECTOR.  ♦ 

***4**4***444***4«***44**4*******44****4***********************4*4 

R12X=X2-X1 

R12Y=Y2-Y1 

R12Z=Z2-21 

DEN0M=SqRT(X2*X2+Y2*Y2+4.0*FM*FM) 

VNX=-X2/DEN0M 

VNY=-Y2/DEH0M 

VNZ=2.0*FM/DEN0M 

TM11=1.0-2.0*VNX*VNX 

TM12=-2.04VNX*VNY 

TM13=-2.0*VNX*VNZ 

TM21=TM12 

TM22= 1.0-2. 0* VNY*VNY 
TM23=-2.0*VNY*VNZ 
TM31=TM13 
TM32=TM23 

TM33=1.0-2.0*VNZ4VNZ 

*****************m******************************»»*****0********** 

*  COMPUTE  THE  COORDINATES  IN  THE  APERTURE  PLANE  OF  THE  MAIN  ♦ 

*  REFLECTOR .  * 

**444i4«44****44*******4********4**4*********4**«***4***4**4i*4****4 

DEH0M=TM31*R12X+TM32*R12Y4-TM33*R12Z 

V1=(TM21*R12X+TM22*R12Y+TM23*R12Z)/DEN0M 

V2=(TM11*R12X+TM12*R12Y+TM13*R12Z)/DEN0M 

Z3=ZM+XE*XE/(4.0*FM) 

Y3=Y2+(Z3-Z2)*V1 

X3=X2+(Z3-Z2)*V2 

RETURN 

END 


1 


SUBROUTINE  PEKPOW(MN,NE,X,Y,XT,YT.ASA,V,ZA,FS,ZS,FM,ZM.XB,XC.XD 
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2 

3  C 

4  C 
s  C 

6  C 

7  C 
6  C 
9  C 

10  C 
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13  C 

14  C 
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16  C 

17  C 
IS  C 

19  C 
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23  C 
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26  C 

29  C 
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31  C 

32  C 

33  C 

34  C 

35  C 

36  C 
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39  C 

40  C 

41 

42 

43 

44 

45 

46  C 

47  C 
46  C 

49  C 

50 


k  .XE.IRC.ISIH.RPEAK.PMAX) 


**4*4**4*4**4*************4*******»4*444*44*4*4***44**4**4«44***« 

THIS  SUBROUTINE  COMPUTES  THE  PEAK  POWER  VALUE  AND  THE  ANGULAR  * 
LOCATION .  • 

«**•**«*««**•****•*««*«*«•«««*«*•«••*•*«*««******«***«*««*««««*«• 

INPUTS :  ♦ 

4 

NN  -  MAXIMUM  ARRAY  DIMENSION.  * 

NE  -  TOTAL  NUMBER  OF  ELEMENTS.  * 

X(MN)  -  COORDINATES  OF  THE  CENTER  OF  EACH  ELEMENT.  * 

Y(MH)  ♦ 

XT(MN)  -  COORDINATES  OF  THE  CENTER  OF  EACH  ELEMENT  PROJECTED  • 
YT(MN)  ONTO  THE  MAIN  REFLECTOR  APERTURE  PLANE.  ♦ 

ASA(MN)  -  ARRAY  OF  CORRECTION  PHASE  ANGLES  IN  DEGREES.  * 

V(MN)  -  ARRAY  OF  TOTAL  PHASE  CORRECTION  INCLUDING  PATH  ♦ 

LENGTH  ERROR  CORRECTION.  ♦ 

ZA  -  ALONG  AXIS  ARRAY  LOCATION  IN  WAVELENGTHS.  * 

FS  -  FOCAL  LENGTH  OF  SUBREFLECTOR  IN  WAVELENGTHS.  ♦ 

ZS  -  ALONG  AXIS  SUBREFLECTOR  LOCATION  IN  WAVELENGTHS.  * 

FM  -  FOCAL  LENGTH  OF  MAIN  REFLECTOR  IN  WAVELENGTHS.  ♦ 

ZM  -  ALONG  AXIS  MAIN  REFLECTOR  LOCATION  IN  WAVELENGTHS.  * 

XB  -  MINIMUM  LATERAL  OFFSET  OF  SUBREFLECTOR  GEOMETRY  IN  * 

WAVELENGTHS .  • 

XC  -  MAXIMUM  LATERAL  OFFSET  OF  SUBREFLECTOR  GEOMETRY  IN  • 

WAVELENGTHS .  » 

XD  -  MINIMUM  LATERAL  OFFSET  OF  MAIN  REFLECTOR  GEOMETRY  IN  ♦ 
WAVELENGTHS .  * 

XE  -  MAXIMUM  LATERAL  OFFSET  OF  MAIN  REFLECTOR  GEOMETRY  IN  • 
WAVELENGTHS .  * 

IRC(MN)  -  INTEGER  ARRAY  USED  TO  MARK  ANY  RAYS  WHICH  SPILL  ♦ 
OVER  THE  SUBREFLECTOR.  ♦ 

ISIN  -  ISIN=0  MEANS  NO  CORRECTION  IS  PERFORMED  AND  ISIN=1  ♦ 

MEANS  FULL  CORRECTION  OF  PHASE  ERROR.  * 

4 

OUTPUTS :  4 

4 

RPEAK  -  PEAK  POWER  ANGULAR  LOCATION  IN  DEGREES.  4 

PMAX  -  PEAK  POWER  LEVEL.  4 


444444444444444444444444444444444444444444444444444444444444444444 

REAL44  X(MN) , Y(MN) .XT(MN) .YT(MN) . ASA(MN) . V(MN) 

INTEGER44  IRC(MN) 

PI=3. 141593 

TP=.628318S3E+01 

RAD=.17453293E-01 

444444444444444444444444444444444444444444444444444444444444444444 

4  TRANSFORM  ALL  ARRAY  COORDINATES  TO  THE  MAIN  REFLECTOR  APERTURE  4 
4  PLANE.  4 

444444444444444444444444444444444444444444444444444444444444444444 

2AP=ZM+(1/(44FM))4XE*XE 
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51 

52 

53 

54 

55 

56 

57 

58 

59 

60  C 

61  C 

62  C 

63  C 

64 

65 

66 

67 

68 

69 

70 

71  C 

72  C 

73  C 

74 

75 

76  C 

77  C 

78  C 

79 

80 
81 
82 
83 
64 

85 

86 

87  C 

88  C 

89  C 

90 

91  C 

92  C 

93  C 

94 

95 

96  C 

97  C 

98  C 

99 


PL0=2*ABS(ZS-ZA)>ABS(ZA-ZM)+ABS(ZAP-ZM) 

XCS=(XB+XC)/2.0 

SRADUS=ABS(XC-XCS) 

DO  1  1=1, BE 
XO=X(I) 

YO=y(I) 

ZO=ZA 

CALL  RAYATS(XO,YO.ZO.ASA(I),FS.ZS.FM.ZM.XE.X1.Y1,Z1 
k  .X2.Y2,Z2,X3,Y3.Z3) 

«**«***«**«**********«*******8*8**** ****** **4* ********»8*«****** 

*  CHECK  TO  SEE  IF  THE  RAY  MISSES  THE  SUBREFLECTOR  OR  THE  MAIH  * 

*  REFLECTOR .  * 

**************************************************************** 

DX=X1-XCS 

RAY=SQRT(DX*DX+Y14Y1) 

IF  (RAY  .LE.  SRADUS)  THEN 
IRC(I)=0 
ELSE 

IRC(I)=1 
END  IF 

**************************************************************** 

*  STORE  THE  APERTURE  COORDINATES.  * 

**************************************************************** 

XT(I)=X3 

YT(I)=Y3 

**************************************************************** 

*  COMPUTE  AND  STORE  THE  APERTURE  SCAN  ANGLE.  * 

**************************************************************** 

T=X3-X2 

B=Z3-Z2 

IF  (T  .EQ.  0.0)  THEN 
ARG=0.0 
ELSE 
ARG=T/B 
END  IF 

ASCAN=ATAN(ARG) 

**************************************************************** 

*  COMPUTE  THE  PATH  LENGTH  ERROR.  * 

**************************************************************** 
CALL  PATHLN(XO.YO,ZO.X1,Y1.Z1.X2,Y2.Z2,X3.Y3,Z3.PL) 
**************************************************************** 

*  GENERATE  THE  PATH  LENGTH  PHASE.  * 

**************************************************************** 
PLP=TP*PL 

PLE=PL-PLO 

**************************************************************** 

*  GENERATE  THE  SCAN  PHASE.  * 

**************************************************************** 

SCPH=-TP*SIN(ASCAN)*XT(I) 
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100  C 

101  C 

102  C 

103 

IM  C 

105  C 

106  C 

107 
103 
109 

no 

111 

112  1 

113  C 

114  C 

115  C 

116 
117 
113  C 

119  C 

120  C 

121 
122 

123 

124 

125 

126 
127 
123 

129 

130 

131 

132 

133 

134 

135  2 

136 

137 


***««*********33**3**«*****3****3******3***** *4**3* ************* 

*  GENERATE  THE  PATH  LENGTH  CORRECTION  PHASE.  « 

**************************************************************** 

PLCP=-TP*PLE 

**************************************************************** 

*  GENERATE  THE  ELEMENT  PHASE.  * 

**************************************************************** 

IF  (ISIN  .EQ.  0)  THEN 
V(I)=SCPH+PLP 
ELSE 

V(I)=SCPH+PLP+PLCP 
END  IF 
CONTINUE 

*******************************************************«********** 

*  SEARCH  RANGE  LIMITS  IN  RADIAMS .  * 

*******************************************»********************** 

A=-RAD*10.0 

B=RAD*17.0 

****************************************************************** 

*  FIND  THE  PEAK  OF  THE  MAIN  BEAM.  * 

****************************************************************** 

HMAX-378 

RsO.O 

P=0.0 

DR=(B-A)/HMAX 
PMAX=-1.0 
RPEAK*0.0 
DO  2  I=0,NMAX-1 
R=A+I*DR 

CALL  AFPAT(MN,NE,XT,YT,V,R.P,PAT) 

IF  (PAT  .GT.  PMAX)  THEN 
PMAX=PAT 
RPEAK=R 
ELSE 
END  IF 
CONTINUE 
RETURN 
END 


1  SUBROUTINE  AFPAT(MB.NE,X,Y,V,T,P,PAT) 

2  C  *****************************************************************; 

3  C  *  THIS  SUBROUTINE  COMPUTES  THE  SQUARE  OF  THE  ARRAY  FACTOR. 

4  C  ***************************************************************** 

5  C  *  INPUTS; 

6  C  * 

7  C  *  MN  -  MAXIMUM  ARRAY  DIMENSION. 

8  C  *  NE  -  TOTAL  NUMBER  OF  ELEMENTS. 

9  C  *  X(MN)  -  COORDINATES  OF  THE  CENTER  OF  EACH  ELEMENT. 

10  C  *  Y(MN) 

11  C  *  V(MH)  -  ARRAY  OF  TOTAL  PHASE  CORRECTION  INCLUDING  PATH 
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12  C  *  LEKGTH  ERROR  CORRECTIOS . 

13  C  *  T  -  THETA  IN  RADIANS. 

14  C  *  P  -  PHI  IH  RADIANS. 

15  C  ♦ 

16  C  •  OUTPUT: 

17  C  * 

18  C  ♦  PAT  -  SQUARE  OF  ARRAY  FACTOR. 

19  C  t**************************************************************** 

20  REAL*4  X(MN),Y(KN) .V(MN).KX,KY 

21  TP=.628318S3E+01 

22  CP=COS(P) 

23  SP=SIN(P) 

24  ST=SIN(T) 

25  KX=TP*ST*CP 

26  KY=TP*ST*SP 

27  SUM1=0.0 

28  SUM2=0.0 

29  DO  1  1=1. NE 

30  ARG=V(I)+KX*X(I)+KY*Y(I) 

31  SUM1=SUM1+C0S(ARG) 

32  SUM2=SUM2+5IN(ARG) 

33  1  CONTINUE 

34  PAT=SUM1*SUM1+SUI12*SUH2 

35  RETURN 

36  END 
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